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Abstract. 

We present optical to radio spectral energy distribution fitting of the host galaxies of four long 
gamma-ray bursts: 980703, 000210, 000418 and 010222, which were detected at submillimetre 
and/or radio wavelengths. We find that only very young starburst galaxy models are consistent with 
the data having both blue optical colors and a pronounced submm emission. For each host we are 
able to construct a model consistent with the short- and long-wavelength parts of the spectra. We 
find galaxy ages ranging from 0.09 to 2.0 Gyrs and star formation rates ranging from 138 to 380 
M yr _1 . 
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INTRODUCTION 

Observations of gamma-ray burst (GRB) host galaxies have attracted a lot of interest 
since the first interpretation of a constant extended source at a GRB position as the host 
galaxy (U, @]. A detailed description of the GRB environment can provide important 
constraints on GRB physics [see[3i for a review] and can be used for statistical studies 
of starburst galaxies. 

The host galaxies of the long GRBs 000210, 000418 and 010222 are somewhat special 
because they are the only hosts, among ~30 targeted by SCUBA [4], showing submil- 
limetre emission. The hosts of GRBs 980703 and 000418 were also firmly detected at ra- 
dio wavelengths [H, S 0] whereas those of GRBs 000210 and 010222 were only weakly 
detected [7]. The host of GRB 010222 was also detected at millimetre wavelengths [Jstl . 
Both submillimetre and radio detections may indicate high star formation rates (SFRs), 
of the order of several hundreds solar masses per year (M yr _1 ), if the whole emission 
is powered by a starburst, namely by dust at submillimetre and by supernova remnants 
at radio wavelengths. Moreover, a spectral energy distribution (SED) fittin g ap plied only 
to optical/near-infrared data agreed reasonably with starburst templates D9j,|lO|,LLl|]. This 
picture is consistent with the hypernova GRB model [j l2l[l3|l . 



However, the starburst scenario of GRB hosts is more complicated because, unlike 
red, dusty, sturbursting submillimetre galaxies [14], GRB hosts have been found to 
exhibit blue optical colors [@, Qii 15, [3. The submillimetre/radio fluxes were also 
underestimated by any SED template fitted to the optical data. 

Mid-infrared Spitzer observations of the hosts of GRBs 980703 and 010222 could 
not bridge the optical and submillimetre data and provide an explanation of the unusual 



TABLE 1. Ages and total SFRs of GRB host galaxies. All four are young starbursts. 
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Defined as the time since the beginning of the galaxy evolution. 

Errors are statistical at 1(7 level calculated assuming that the template fits to the data. 
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SED behaviour and a clue on dust properties because only the former was detected in 
the bluest 4.5 /im Spitzer passband fl 1 VI . 1 1 811 . 

In this contribution we attempt to find reasonable full SED fits for those four GRB 
host galaxies and constrain their ages and SFRs. This approach seems to be a promising 
tool to constrain the dust properties and the obscured star formation rates as shown, for 
example, by Priddey et al. Ill9n . We use a cosmological model characterized by Hq = 70 
km s" 1 Mpc" 1 , Q A = 0.7 and Q M = 0.3. 



SED MODELLING 

In order to model the SEDs of GRB hosts we used the GRASIL0 software described by 
Silva et al. Ir20ll . It is a numerical code that calculates the spectrum of a galaxy by means 
of a radiative transfer method, applied to photons produced by a stellar population, and 
reprocessed by dust. The importance of this model is the fact that it is self-consistent in 
that it fulfills the principle of energy conservation between the energy absorbed by dust 
in the UV/optical wavelengths and the energy re-emitted in the infrared. 

Figured] shows that submm and radio emissions of the GRB hosts discussed here are 
underestimated even for a template corresponding to the nearby ultraluminous infrared 
galaxy (ULIRG) Arp 220 |2(J, which is also too red in the optical compared to GRB 
hosts. This indicates the need for new empirical SED models, which are presented here. 

We performed SED modelling investigating a wide range of GRASIL parameters (see 
Michalowski [|2l|l for a description) We scaled the SED templates to the datapoints and 
chose those which resulted in acceptably small % 2 . Then we calculated the total SFRs 
using the infrared luminosity, integrating between 8 and 1000 /im (rest frame), and 
applying Kennicutt I& SFR(M yr" 1 ) = 4.5 • 10" 44 L 8 _i 00 o(erg s" 1 ). 



|http : / /web . pd. astro . it /granato/grasil/grasil . html] 




FIGURE 1. SEDs of GRB hosts. Solid lines: the young starburst galaxy models calculated using 
GRASIL and consistent with the data. Dotted lines: Arp 220 model [from 20]. Squares: detections with 
errors, in most cases, smaller than the size of the symbols. Arrows: 3(7 upper limits (values marked at the 
base). Hashed columns mark the wavelengths corresponding to optical/near infrared filters, mid-infrared 
Spitzer filters, SCUBA submillimetre bands and the radio domain. Data are from: optical/near-infrared: 
Gorosabel et al. ll9l [Toll . Vreeswijk et al. Il5ll . Galama et al. lfl6ll . Sokolov et al. j2m : mid-infrared: Le 
Floc'h et al. lfl7ll . Castro Ceron et al. Il8ll : submillimetre: our reduction of archival SCUBA data, Tanvir 
et al. [4]; radio: Berger et al. SH, Frail et al. H, Sagar et al. Q 



RESULTS AND DISCUSSION 



For each host we were able to construct an SED that fitted the data reasonably well. The 
corresponding ages of the galaxies and the SFRs are presented in Table Q] and fits are 
shown in Figure [TJ 

All the galaxies discussed are young (ages less than 2 Gyrs). Here, the age of a galaxy 
was defined as the time since the beginning of its evolution, when the stellar population 
starts to build up. The bigger the difference between optical and submillimetre fluxes, 
the lower the needed age was. This is because a younger galaxy has more stars still 
embedded in molecular clouds, so optical light is weaker and dust emission stronger. 

The construction of the SEDs assuming low ages for the galaxies was therefore the 
way to take into account both significant submillimetre/radio emission and the blue 
optical colors of GRB hosts. In young galaxies there are, on one hand, lots of young, hot, 
blue stars, because they have not finished their lives yet. Hence the total optical spectrum 
of the galaxy is blue. On the other hand, the majority of the stars still resides in dense 
molecular clouds, so a significant part of the energy is absorbed and re-emitted. This 
increases the dust emission. GRBs reside in molecular clouds and it was found that gas 
column densities derived from X-ray afterglows in a sample of 8 GRBs (including GRB 
980703, discussed here) were in the range corresponding to the column densities of giant 
molecular clouds in the Milky Way OQTI . A similar conclusion for high-redshift GRBs 
was recently drawn by Jakobsson et al. ||3~lll by means of modelling Lya absorption 
features. 

Based on SEDs fitted to optical data, Gorosabel et al. |@, and Christensen et al. 
lUlll derived the ages of the starbursts in the hosts of GRBs 980703, 000210, 000418. 
Our estimates are at most a factor of two larger. Given that we calculated the time from 
the beginning of the galaxy evolution, not the beginning of the starburst, our values agree 
with those derived by these authors. Sokolov et al. 112 811 derived ages of both old stellar 
populations and starbursts. Our estimation for GRB 980703 agrees with the age of the 
old component, which is conceptually closer to our definition of the galaxy age. 

Table Q] also presents the SFRs derived from the SED fits. They are all high, of 
the order of several hundreds solar masses per year. Statistical la errors, calculated 
assuming that the model represents the data, are also given in Tabled! On top of that, 
one should add the 30% uncertainty of the relation between infrared luminosities and 
SFRs 

In TableCQour SFR estimations are compared with those derived using three other SFR 
indicators: mid-infrared, submillimetre and radio emission. Our results are consistent 
with radio-derived SFRs |7D. This is because the calibration of SFR to radio flux requires 
choosing only two parameters [a normalization factor and a spectral index, see l32ll. 
which are relatively well constrained. Our values also agree with the upper limits derived 
by Castro Ceron et al. II 1 311 using the template of Arp 220. This approach appears to 
be more self-consistent and reliable than analysing only 24 /im datapoints, which is 
strongly dependent on the unknown shape of the infrared part of SED. For example, 
this method applied by Le Floc'h et al. [17] gave upper limits inconsistent with our 
results. Finally, Berger et al. [7] obtained higher SFRs based on submm alone. 

The mid-infrared (5 }im < X < 40 /im) shape of the SEDs was not constrained due to 
a lack of targeted observations or detections. Hence, all SEDs shown on Figure CD have 



different mid-infrared behaviour depending on type, mixture and sizes of dust grains 
used. Any data in this region would be useful to learn more about the dust properties in 
these galaxies. 

This work provides strong evidence that some GRB hosts are young starbursts. We 
have proposed an explanation why they have both blue optical colors and significant sub- 
millimetre emission unlike red, dusty submillimetre galaxies UM- The SED approach is 
self-consistent and makes use of all the data available. The research presented here will 
be continued to build up a larger galaxy sample analysed in this way. Other properties of 
the host galaxies, the robustness of the technique and a more accurate treatment of the 
errors are discussed in Michalowski et al. ||33l. l34fl . 
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